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ABSTRACT 

The Lower Cretaceous Mannville Group in west-central Saskatchewan 
jis wedge-shaped, thickening towards the north. It lies unconformably on 
Devonian strata in the north and in the south mainly on Mississippian 
beds which form an east-west trending ridge. 

The basal Dina Member of the Mannville is thickest in erosional 
valleys on the sub-Cretaceous topographic surface and thins over ero- 
sional highs. Following deposition of the Dina sandstone, relief on the 
depositional surface was subdued. The Cummings Member, which has a thin 
shale unit at the base and peeore: sandy towards the top, was deposited 
by the Clearwater Sea over this surface. The Dina and the Cummings Mem- 
bers do not extend farther south than Township 40. 

A coal seam ("seam 3") was Henecived near the top of the Cummings 
Member, indicating a temporary regression of the Clearwater Sea. An- 
other major coal seam ("seam 8") was deposited during late Mannville 
time. This latter seam is absent, however, between Townships 37 and 39, 
probably due to non-deposition. 

The Battleford Arch passes through the central part of the tnesis 
‘area. It is present as a topographic high which probably formed a depo- 
sitional barrier during deposition of the Dina and the Cummings Members. 
The absence of seam 8 in Townships 37 to 39 is probably also due to this 
high. 

The coal in the Mannville Group is of no commercial value at the 
present time. It is scattered throughout the Lower Cretaceous sediments. 
Resources in place are estimated to be approximately 19,4060 million 


(short) tons (approximately 17,600 million metric tons). 
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Chapter | 
INTRODUCTION 


A. LOCATION OF STUDY AREA AND SOURCES OF INFORMATION 

The area studied in this thesis occupies approximately 8,700 square 
miles (22,560 square kilometers) in west-central Saskatchewan, extending 
from Range 18 west of the Third Meridian on the east to the Alberta- 
Saskatchewan border on the west, with the southern and northern bound- 
aries including Townships 29 and 50, respectively (figure 1.1). Data 
were obtained from geophysical logs of 236 wells (figure 1.2) with 115 
wells having coal occurrences. Another 101 wells were examined ina 
small area (see figure 8.3) to obtain closely-spaced data for the vario- 


gram study. 


B. OIL AND GAS FIELDS AND COAL IN THE STUDY AREA 

The major oi] fields in the area are: Dodsland, Coleville-Smiley, 
North Hoosier and Lloydminster. The productive zone of the former two 
fields is the Viking Formation (Lower Cretaceous). In the North Hoosier 
field iqndersley area) oil is found in the Middle Bakken Sandstone 
(Mississippian), Deville Member (Detrital Zone) and Dina Member (Basal 
Blairmore Sand), of Early Cretaceous age (White, 1974). The Residual 
Zone Pool, as it is called, in the vicinity of Section 11-31-28 west 
of the Third Meridian, is the only producing reservoir in the Deville 
Formation. Oil is commercially produced in the North Hoosier Basal 
Blairmore Sand Pool in Sections 10 to 15 and 22 to 24, Township 32, 


Range 27 west of the Third Meridian. The reservoir is several square 
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Fig 14.1 LOCATION MAP. 
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miles in areal extent and is capped by an impervious layer of shale or 
mudstone and underlain by coal (seam 6). Both pools produce heavy oi] 
with gravity between 13° and 15° API (specific gravity, 0.979 to 0.966). 
Reasoner and Hunt (1954) suggested that tere source for the Bakken and 
Mannville accumulation was the Mannville sediments. 

| The North Hoosier field produces 28 percent of Saskatchewan's heavy 
oil. The remaining 72 percent of heavy oi] production is derived from 
sands of the Mannville Group {mainly the Sparky sand) in the vicinity of 
Lloydminster. 

Altogether, ten coal seams are present within the Mannville Group, 
but only two (seams 3 and 8) have regional significance. Seam 3 was 
deposited near the top of the Cummings Member and extends as far south 
as the member itself (Township 40). Seam 8 covers the northern and 
southern part of the study area and was deposited on top of the Sparky 
sand (equivalent to the upper part of the Borradaile Member). The re- 
maining eight seams have very limited local extent; two of which (seams 


1 and 2) were encountered only once in the 236 wells studied. 


neC fa-PURPOSECOR: THE THESIS 
Maycock (1967) discussed the geology of the Lower Cretaceous sedi- 
ments in west-central Saskatchewan, and between Townships 29 and 40, he 
referred to the Mannville sediments as "Undifferentiated Mannville 
Group". On the other hand, Nauss (1945) subdivided the Mannville into 
six members in the Vermilion area, approximately 36 miles west of the 


“northwest corner of the present study area. One of the purposes of this 


thesis is to find the boundary between these two classifications and what 
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constitutes this boundary. Correlation between the two classifications 
based on lithology is evidently inapplicable. Because there is a lack 
Of f0ss 1iS,ouSiIng fossiis for correlation is not possible at the present 
time. In this thesis, the writer had chosen coal as a means to explain 
some of the stratigraphic problems. A number of coal seams do occur 
both in the northerm and in the southern part of the study area, provid- 
ing a good tool for correlation purposes. The writer had also used coal 
to gain some insight into the depositional environments. Finally, in 
order to give an idea of how much coal is present in the Mannville sedi- 
ments studied, a resources estimate was carried out. To summarize, this 


thesis project was undertaken in order to: 


(1) determine the areal relationships between Mayceck's and Nauss' 
Stratigraphic subdivisions, 

(2) attempt to determine the causes of the facies changes between 
Maycock's and Nauss' areas, 

(3) study the distribution of coal in the Mannville Group of western 
| Saskatchewan and to determine whether coal distribution can aid in re- 
constructing depositional environments in the Mannville Group, and 


(4) estimate the in-place coal resources of the area mapped. 
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Cnapter 2 
PREVIOUS WORK 


Because of the hydrocarbon potential of Lower Cretaceous strata 
in the Lloydminster area, the Mannville Group has drawn much attention 
recently. Earlier work which deals with the Mannville Group in the 
Study area is limited srther in areal extent or in geological aspect. 
Reasoner and Hunt (1954) discussed the geology of the Coleville-Smiley 
oil field; Price (1955) extended the term Mannville from Alberta into 
Saskatchewan. More recently, Maycock (1967) studied the area between 
Townships 18 and 40, Ranges 15 to 29 west of the Third Meridian in some 
detail. Studies on local oil and gas pools have also been carried out 
by Wickenden (1948) and White (1969, 1972, 1974). Other work that deals 
with neighbouring areas include Nauss' (1945) subdivision of the Lower 
Cretaceous beds in the Vermilion area based on the presence or absence 
of dark minerals, Badgley's (1952) description of the Mannville Group 
in central Alberta, Glaister's (1959) correlation of the Lower Creta- 
ceous of southern Alberta and adjoining areas, Williams’ (1963) reclas~ 
sification of the Mannville Group in central Alberta based on petro- 
drapheel eee raph IE ey I Uae and Mellon's (1967) petrographic study 
of the Mannville Group in Alberta. The most recent work is that of 
Christopher (1974) in southwestern Saskatchewan in which he recognized 
a new formation (Success Formation) at the base of the group, below the 
Cantuar and Pense Formations used by Price (1955) and Maycock (1967). 

The name Mannville was first used by Nauss (1945) for a predomi- 


nantly non-marine formation in the Vermilion area, lying between the 
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sub-Cretaceous unconformity and the marine Joli Fou shale. The term has 
been used for the Lower Cretaceous sediments within the Alberta Plains 
for some time, yet correlative beds are still popularly referred to as 
Blairmore within the Saskatchewan subsurface. These rocks, however, are 
both Heographi ca uly and lithologically different from those of the type 
Blairmore of the Alberta Foothills. Price (1955) introduced the term 
Mannville into Saskatchewan and subsequently (1963) showed the areal ex- 
tent of these strata within the Interior Plains (figure 2.1). 

The first estimate of coal resources in Canada was prepared by 
Dowling (1913) for the 12th International Geological Conarese: His es- 
timate of Canada's pezolrtces was 1.3 trillion (short) tons, based on coal 
seams with a minimum thickness of one foot to a maximum depth of four 
thousand feet. Mackay (1947) reported 99 billion tons of coal in-place 
in his study for the Royal Commission on Coal in 1946. The reason for 
Mackay's lower ene is that he used seam thickness and depth of over- 
burden that varied from coal field to coal field and did not include 
seams less than three feet thick. Latour and Chrismas (1970) made a 
more detailed survey on western Canadian coal, using tne following pa- 
rameters: | 

(1) seams less than three feet thick and rock partings that are one 
foot thick or thicker are not included in calculating the tonnage, 

(2) only coal seams that are within one hundred and fifty feet from 
the surface are considered. 

They estimated a total of 118.7 billion tons of coal in-place in 
western Canada. Of this, 9.8 billion tons (8.3%) are considered measured 


resources, 50.3 billion tons (42.4%) indicated and 58.6 billion tons 
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(49.3%) inferred. Among the three western provinces, Saskatchewan has 
10.1% of the total coal resources in Late Cretaceous to Tertiary (Raven- 
scrag Formation) beds in the south of the province. 

Yurko (1976) surveyed the deep Cretaceous coal resources in the 
Alberta Plains and estimated a total of 628 billion tons for seams two 
feet thick or thicker in the Mannville Group. This represents 60% of 
the estimated total deep coal resources in the Cretaceous System in 
Alberta. The method used was similar to that used by the writer for 
this study. 

Other than these regional estimates for the western provinces, no 
data has been published on the estimate of the coal resources in the 
study area. Coal seams in the Mannville Group are generaliy thin and 
about 1,700 to 2,000 feet below the surface. This coal is of no commer- 
cial significance at the present time, but could become important if 


techniques for underground "in situ" gasification are developed. 


-iivnttet alee! , 25 yet Wisteew ae 


noweti? wantdant a: Peertr yy B: nha z yas " 
ge MONG: git 1% inion point Gk 
Od, Pi wOTAE RS Tein hee ora qgeed” ae aay 2 uret 
id sglhee totuchos Nate (Rs snd a basanites hae ene St 
Lua gihegsagyy & rea ivan S(T? NN pas I mast * <1 
"misren Sew. Maasai: mt eile egis ellen “ehh 
. Sob-neri vr ond ed Weenies. a “WAR! ee wth a fs 
ROM TD ec Men saicetoys sealagn nah nity antlte 
* oft ~ Tifa? Seeger Jap 5) ate a nan en ant he ‘gatas boa it 
~ “ting Ade” y gene’ Gp eae a gt ‘amiga $00) aye 2 ite, 
-cameena- 4a: eo Bia teow agli? - senthos sits wape sh MG betiaial ia ) 
*s “Ushadveelint. pene? et aew, bat Sat Serna ae 
benotovah 44 ee hatte wit ison 


a 


: 
wf 
‘ A " 7 
‘ t = 
« : 
4 Re 
, . 
“s 
i] ' 
: > ; - cal 
7 a ~~ . -~ 
' vi 
“ - \ 
7 ei f > p > _ 
- - a 
; s 
be a 
aM A = 
ts a2 a 
‘ 
* ' es 
‘ 5 5 
= > -~ «- > 
“; Pe é 
s 
a 
= ey Oh: ; . Soe 7 
’ rf 
~ oe. a 7 
i -_ 
- * ta - al ~ 


Chapter 3 
METHOD OF STUDY 


A. AVAILABILITY OF INFORMATION 

Data for this study were obtained from geophysical borehole logs, 
which were run primarily for oi] exploration. Therefore, a complete 
suite of logs favourable for coal identification was not always avail- 
able for each borehole. All 236 wells used in the regional study have 
at least a resistivity log; 95 of the 236 wells have a sonic log; 25 
have a neutron log; 5 have a density log; 7 have both sonic and density 
logs; 6 have both sonic and neutron logs, and | has density and neutron 
logs. In the remaining 97 wells, data were derived solely from the re- 
sistivity log. The resistivity log always includes a spontaneous po- 
feieial (SP) log, and the sonic, density and neutron each include a 


gamma log. 


B. LOG CHARACTERISTICS OF THE MANNVILLE GROUP 

The top of the Mannville Group is indicated by the first increase 
in resistivity on the electrical Jog below the Viking sandstone. This 
increase is always accompanied by a corresponding negative defiection 
of the spontaneous potential (SP) curve (figure 3.1). North of Township 
40, approximately 160 feet below the top of the Mannville, a high resis- 
tivity zone with high negative SP defines the Sparky sandstone (equiva- 
lent to the upper part of the Borradaile Member). The thickness of this 
sandstone varies and the base is not always distinct. Therefore, only 
the top was picked for correlation purposes. 


The Cummings and Dina constitute the basal members of the Mannville 
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LOG CHARACTERISTICS of the MANNVILLE GROUP. 


Fig 3.1 
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Group. Within the thesis area, the Cummings Member is a sandstone with 
a high negative SP value and very low resistivity (figure 3.1). These 
characteristics also apply to the Dina Member. Separating the two is 

a positive deflection in the SP curve at the base of the Cummings. 

The Sparky sandstone and the Cummings and Dina Members lose their 
distinctiveness from north to south, and south of Township 40 they dis- 
appear completely (figure 3.2). South of Township 40, no attempt was 
made to differentiate the Mannville sediments into individual members, 
and the writer has adopted Maycock's (1967) designation as "Undifferen- 
tiated Mannville Group". 

The base of Mannville is easier to pick in the north than in the 
south. North of Township 35, the Mannville is underlain unconformably 
by Devonian carbonate rocks with very high resistivity. Since the Dina 
Member has a very low resistivity, the sharp break in the resistivity 
curve indicates the position of the sub-Cretaceous unconformity (figure 
3.1). South of Township 35, Lower Cretaceous rocks are underlain in 
places by Mississippian beds, and the boundary between Mississippian and 
Cretaceous formations is not always clear on the geophysical logs. The 
Deville Member, if present, has a high negative SP signature, but unlike 
the Dina, it also has a high resistivity. Mississippian limestones, 
like those of the Devonian, always show high resistivity and high nega- 
tive SP, and it is difficult to distinguish between them and the Deville 
Member. The Mississippian Bakken Formation is subdivided into three 
(upper, middle and lower) members, each characterized by a distinctive 
"spike" on the resistivity log. The upper and lower members consist of 
shale, separated by a sandstone which in places is oil-bearing (figure 
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Coal has a very characteristic signature on geophysical logs and 
a detailed discussion on this topic has been published by Bond et al. 
(1966. (see also Yurko, 1976). The following account is a summary of 
some of the diagnostic properties of coal in borehole interpretation 
(see also figure 3.3). Figure 3.4 shows a composite log of the relative 
position of the coal seams in the study area. 

To Density. 

Coal densities range from 1.29 to 1.8 gm/cc, dependine on the rank, 
depth of burial, and ash and moisture content. Coal in the Mannville 
Group is approximately 1,500 feet underground and has a density of 1.8 
gm/cc or less. Therefore, when a density log is available, coal is 
easily recognized by a very low density response as compared to adjacent 
sandstones and shales. 

2 SOIC Lod 

The sonic log records the interval transit time (aT) in microseconds 
per foot, the time required for a compressional sound wave to travel 
across one foot of formation. The AT value for coal is usuaily greater 
than that for shale or sandstone and varies inversely with degree of 
compaction. That is, coal seams closer to the surface have a longer in- 
terval transit time than deeper ones. Tertiary lignite coals of south- 
western Alberta and Saskatchewan have AT values between 140-170 micro- 


seconds per foot, and coals in the Mannville (this study) have AT values 


between 120 and 160 microseconds per foot. 
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Fig 3.4 COMPOSITE LOG SHOWING RELATIVE POSITION of SEAMS. 
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34 Neutron-Log 


The neutron log primarily measures the effectiveness cf a formation 
to moderate or slow down a beam of fast neutrons. Essentially, the mod- 
eration of the fast neutrons is caused by hydrogen atoms in pore fluids; 
however, carbon is also an effective moderator and consequently coal is 
indicated by low readings on the neutron curve. Unlike sonic or density 
logs, the neutron log does not have a definite value range to indicate 
the presence of coal, therefore in using the neutron log to detect coal 
seams, one has to look for relative low values in the curve rather than 
absolute values of the neutron count. 

4. Natural Gamma Log 

Geophysical logs record only the natural gamma radiation in a well 
bore, which is produced by potassium, uranium and thorium in sedimentary 
rocks (mainly in shale}. Coal, as an organic sediment, contains very 
small quantities of these elements, therefore has very low radioactivity. 
Dee UReSISbiv Icy og 

A relatively high, electrical resistivity is characteristic of clean 
coal; however, the actual value recorded on a log depends on the ash and 
moisture content of the coal. Dirty coal, carbonaceous material and even 
woody fragments may produce low resistivity values on the electrical log. 

Coal beds are indicated by high resistivity spikes on the lateral 
curve with adjacent very low readings in the blind zones immediately be- 
low the coal. On normal resistivity curves, if the coal beds are thinner 
than the electrode spacing, reversals (low readings) generally occur op- 


posite coal beds. 
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No definite response is registered opposite a coal bed although in 
the Mannville Group most seams show a slightly negative SP. The spon- 
taneous potential registered on the SP curve depends mainly on the re- 
sistivity and thickness of the coal bed, the resistivity and the diameter 


of the zone contaminated by mud filtrate. 


D. RELIABILITY INDEX 
With a combination of resistivity and other logs, coal can be iden- 
tified with a degree of certainty which can be indicated By a numeric 
reliability index (see Yurko, 1976, pp. 11-12). In this thesis, the 
writer chose to use reliability indexes 1, 3 and 5, corresponding to 
indexes 1, 3 and 5 used by Yurko (1976), with the following criteria 
used in assigning these numbers. 
(1) Reliability index 1: coal picked from an electrical log and 
at least two of the following specific curves: 
(a) natural gamma, 
(b) sonic (acoustic), 
(c) density, 
(d) neutron. 
Probability of coal actually being present is estimated to be greater 
than 902. 
(2) Reliability index 3: coal picked from an electric log, and 
neutron or natura}? gamma log. As mentioned previously, no specific 
| value range is assigned to the neutron count for coal beds and picks 
depend solely on the shape of the curve. The probability of coal actu- 


ally being present is estimated to be between 50 and 75%. 
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(3) Reliability index 5: coal picked from an electrical log only. 
In this case, the probability of coal actually being present drops to 
25% or less. 
| The reasons for not using reliability indexes 2 and 4 as defined 
by Yurko (1976) are that the geophysical logs used in this study are of 
relatively good quality and the log characteristics of coal are rela- 
tively well defined. Using the tops of the Sparky, Cummings and Dina 
Members as markers, coal seams within the Mannville Group were correlated 


and each seam numbered from 1 to 10 in ascending order of the seams. 


E. GENERAL DISCUSSION OF COMPUTER APPLICATIONS 

The data obtained from the geophysical logs were put into a computer 
file from which subsequent data files were built to meet the needs of. 
particular application programs. The plotting of contour maps was done 
by using the "CGRIDI" and "CONTUR" plotting subroutines available from 
the Computing Services program library at the University of Alberta. The 
CGRID1] routine accepts as input the irregularly distributed borehole data 
and produces a regular data grid. The number of grid points (i.e., the 
size of the grid cell) over a map area is controlled by the user; to ob- 
tain good results it is advisable to have at least one data point per 
grid cell. The output from CGRID] is used as the input to CONTUR which 
prepares the maps for subsequent plotting on the Calcomp Plotter. CONTUR 
allows the user to control the number of contours, contour interval and 
the labelling of the contours. The title, scale and axis of the maps 


were done using other plotting subroutines, also available from the sub- 


routine library. 
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Well location maps were generated by using the Generalized Map Post- 
ing Routine (GPR) from the Western Canada Coal Resources Data Base Pro- 
ject. 

A small area seven miles east of Lloydminster was chosen by the 
writer for a detailed study of areal variance of coal seam thickness. 

One hundred and one well logs were studied and the data were used to 
generate statistical tables and variograms in graphical form (see Chap- 
ter 8). A concise account of the theory and application of areal vari- 
ance studies is given by David (1974) and Matheron (1963). A brief 


summary will be presented in Chapter 8. 
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Chapter 4 
STRATIGRAPHY 


A. INTRODUCTION 

Studies by several workers lead to the conclusion that the Mannville 
sediments in the northern plains of Canada can be subdivided into three 
major portions - a marine shaly unit (equivalent to the Clearwater Forma- 
tion) which is sandwiched between a lower relatively homogeneous sandy 
unit (the McMurray Formation) and an upper, more heterogeneous sandy unit 
(the Grand Rapids Formation). 

The lower sandy unit fills in the topographic lows on the sub-Creta- 
ceous erosion surface. In the Interior Plains, the roundness of the 
quartz grains and the absence of dark minerals in this unit suggests a 
high degree of maturity and an eastern source area. Tne overlying shaly 
unit was deposited during the invasion of the boreal Clearwater Sea. 

One of the intentions of the writer is to determine the location of the 
shoreline within the study area. With the withdrawal of the Clearwater 
Sea and the increase in tectonic/igneous activity in British Columbia, 
the upper heterogeneous sandy unit was deposited with its main source 
from the west. As might be expected, this sandy unit has a high content 
of feldspar and volcanic rock fragments. 

Where Mississippian rocks crop out on the sub-Cretaceous surface, 
their weathered products constitute a distinct cherty unit (the Deville 
Member) below the relatively homogeneous sandy McMurray Formation. 

Nauss (1945) first applied the term "Mannville" to Lower Cretaceous 


sediments in the Vermilion area. He subdivided the Mannville Formation 
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into six members as follows in order of increasing age: O'Sullivan, 
Borradaile, Tovell, Islay, Cummings, and Dina (figure 4.1). This sub- 
division can also be applied to the northern plains of Saskatchewan. 

In central Saskatchewan, Maycock (1967) called the correlative sediments 
"Undifferentiated Mannville Group", and further south, these beds inter- 
finger with the Cantuar and the Pense Formations of Price (1963). The 
dividing line between the six-member Mannville Formation as defined by 
Nauss (1945) and the Undifferentiated Mannville Group of Maycock (1967) 
lies within the present study area, in the vicinity of Township 40 (fig- 
ures 3.2 and 4.2). Therefore, the study area can be divided into north- 
ern and southern portions with Township 40 as the boundary for discussion 


purposes. 


_B. THE NORTHERN AREA 

North of Township 40, the Mannville Group within the study area can 
be subdivided into the six members recognized by Nauss (1945) (figures 
4.1 and 4.2). The Dina Member appears as a uniform clean sandstone on 
the SP log (figure 3.1). It fills in the palaeo-topographic low areas, 
~ and pinches out southward towards Township 40 (figure 4.3). After depo- 
sition of the Dina Member, the boreal Clearwater sea transgressed south- 
ward and deposited the Cummings Member which also thins towards the 
south and pinches out in the vicinity of Township 40 (figure 4.4). The 
maximum thickness of the Clearwater occurs in both the northeastern and 
northwestern corners of the map area. The shoreline at that time prob- 
| ably trended east-west through the study area, and to the west turned 


north towards Wainwright where the Cummings is also absent (Nauss, 1945). 
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A coal seam (seam 3) near the top of the Cummings Member is limited to 
Townships 40 to 45, Ranges 24 to 28 west of the Third Meridian (figure 
4.5). The presence of this seam and the sandy lithology of the Cum- 
mings in the study area indicate that the area was probably near shore 
and that a temporary regression of the Clearwater sea occurred during 
the deposition of the member. 

The remaining portion of the Mannville sediments can be subdivided 
into the Islay, Tovell, Borradaile, and O'Sullivan Members (figure 3.2). 
All of them have a southern depositional limit in the vicinity of Town- 
ship 40. A major coal seam (seam 8) at the top of the Borradaile Member 
("Sparky sandstone" in oi] industry terminology) was used as a marker 
horizon. Its distribution is shown in figure 4.6. 

In east-central Alberta, the Dina Member as described by Nauss 
(1945) consists mainly of quartz sandstone containing rounded, frosted 
grains, interbedded with siltstone and shale. Dark minerals are rare. 
The overlying Cummings Member is principally dark grey to black shale. 
As the Clearwater sea regressed, the Islay Member was deposited along 
its beaches (Nauss, 1945, p. 1627). The overlying Tovell, Borradaile 
and O'Sullivan Members have lithologies indicating an alternation of 
deltaic and near shore environments. 

The upper portion of the Borradaile Member is the productive sand 
in the Borradaile oil field. Wilson (1947) noted a regionally widespread 
coal seam at the top of his "Upper Wainwright" unit (equivalent to the 
Borradaile Member), about 130 to 145 feet below the top of the Mannville 
and called it the "Wainwright coal". This same coal seam was picked as 


seam 8 by the writer. 
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In central Alberta, the Mannville Group was divided, in ascending 
order, into the McMurray, Clearwater and Grand Rapids Formations .by 
Badgley (1952). The McMurray Formation, correlative to the Dina Mem- 
ber, is composed of the Deville, Ellerslie and Calcareous Members. The 
Devilie Member is the weathered product of the underlying Mississippian 

-and Devonian rocks. The main facies of the McMurray Formation is the 
Ellerslie Member ("Basal Quartz") which is composed mainly of quartz 
sandstone with minor siltstone, and silty, often carbonaceous shale. 

[t is texturally and compositionally mature, and the sediments have 
gone through several depositional cycles (Williams, 1963). Williams 
et al. (1962, p. 324), using K-Ar dating on the detrital feldspar from 
the McMurray Formation, concluded that the "... K-Ar dates ... indicate 
Precambrian sources on the Canadian Shield ...". 

The top unit of the McMurray Formation is the Calcareous Member 
(Glaister, 1959), characterized by abundant invertebrate remains in- 
cluding pelecypods, gastropods and ostracodes of non-marine and near- 
marine origin (Mellon, 1967). As a lithologic unit, it is also common- 
ly known as the "Ostracode Zone", and has been dated as late Aptian or 
early Albian in age (Glaister, 1959). The Calcareous Member is absent 
in the study area and the Dina Member is overlain by the Cummings Mem-~- 
ber. 

The Clearwater sea transgressed into central Alberta where the 
Clearwater Formation (Williams, 1963) was deposited. The Clearwater- 
McMurray contact, as well as the Clearwater-Grand Rapids contact, are 

| probably diachronous (see Williams, 1963, figure 8; and Mellon, 1967, 


p. 64). For this reason and the fact that the Clearwater Formation 
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cannot be distinguished as a separate entity over much of central Alberta, 
Mellon (1967) combined the Clearwater and the overlying Grand Rapids For- 
mation into the Fort Augustus Formation. 

The Clearwater Formation in the Edmonton area is composed of the 
basal Wabiskaw Member (commonly known as "Glauconitic Sandstone") and 
the overlying dark grey shale. Mellon (1955) suggested that the glauco- 
nitic sandstone is an offshore bar deposit (figure 4.7). Both the Cal- 
careous Member and the "Glauconitic Sandstone" are absent in east-central 
Alberta and in the study area. Nauss (1945, p. 1627) hypothesized that: 
"After the Cummings epoch the sea receded and on its beaches originated 
the glistening sands of the Islay Member." This indicates that the . 
Vermilion area was emerging while deposition of the Clearwater Sea con- 
tinued in central Alberta. If this is the case, Badgley's (1951) cor- 
relation of the "Glauconitic Sandstone" with the Islay Member is ques- 
tionable. 

The Grand Rapids Formation is the uppermost unit of the Mannville 
Group in central Alberta and is mainly of non-marine origin. In the 
Edmonton area, its base interfingers with the Clearwater Shale and is 
diachronous, being older towards the south (Wiliiams, 1963). Williams 
et al. (1962) dated two samples of detrital feldspar from the Grand 
Rapids Formation in Alberta at 120 to 150 and 160 to 220 million years, 
dates which fall within the range of ages from early Mesozoic intrusions 
in the central Cordillera. 

In southern and central Alberta, the Grand Rapids Formation is 
| included in the Fort Augustus Formation by Mellon (1967) and in the 


Upper Mannville by Glaister (1959). Both reported the formation as 
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being mainly of continental origin. According to Mellon (1967, p. 65), 
"Foraminifera are absent from the Brook Stanmore No. 1 well [Sec. 22-30- 
11 W4M], where the entire formation appears to be of non-marine origin, 
and also from the West Viking No. 1 well [Sec. 11-48-15 W4M], although 
the shale at the base of the formation there appears to be marine ...". 
Farther north in the Edmonton area (Twp. 51 to 56) the same writer ob- 
served foraminifera in the lower and central parts of the formation, 
which probably occupy a stratigraphic position similar to that of the 
Clearwater Formation of Williams (1963). The top 140 feet of the for- 
mation include coal-bearing strata interfingered with marine foramin- 
jfera-bearing shales, indicating a fluctuating shoreline. Stili farther 
north (north of Twp. 72), the Fort Augustus Formation has abundant for- 


aminifera in most of the shaly strata and is probably mainly marine in 


origin (Mellon, 1967). 


C. THE SOUTHERN AREA 

The three-fold subdivision of the Mannville Group as defined by 
Williams (1963) or the six-member Mannville Formation of Nauss (1945) 
cannot be extended into the southern portion of the study area. The 
Lower Mannville and the lower part of the Upper Mannville disappear 
eastward into the area where the Cantuar Formation (Price, 1963; May- 
cock, 1967) is recognized. This pinching out results in part from the 
fact that the pre-Mannville surface rises towards the east. 

The Cantuar Formation was established by Price (1963) for the 
| basal Lower Cretaceous sediments in southern Saskatchewan, and the 


Pense Formation for the strata between the Cantuar Formation and the 
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overlying Joli Fou shale. The Pense and the Cantuar together form the 
"basal Cretaceous sandy group" of Price (1963, p. 5). South of the 
Mississippian erosional edge (i.e., south of Twp. 35) in the study area, 
Maycock (1967, p. 13) used the term "Continental Facies" in lieu of 
Cantuar Formation and north of the erosional edge, correlative sediments 
were referred to as "Undifferentiated Mannville Group". Since the Pense 
Formation cannot be recognized within the study area (see also Maycock, 
1967, figure 12) the writer has adopted Maycock's terminology, "Undif- 
ferentiated Mannville Group", to describe the Mannville beds in the 
southern part of the study area. The Undifferentiated Mannville Group 
and the Continental Facies have very similar lithology (Maycock, 1967, 
p. 20) and the following description is applicable to both. 

Maycock (1967) described his Continental Facies as "extremely het- 
erogeneous" and consisting of lithic sandstones associated with argil- 
laceous siltstone, sandstone of variable clay content, lignite, carbo- 
naceous shales and mudstone, together with minor rocks rich in carbon 
filaments. The lithic sandstone he described exhibits calcitic replace- 
ment of the matrix and has relatively simple mineralogy, with the main 
components, in order of importance, being quartz, feldspar, chlorite 
and/or biotite and muscovite. The presence of the less stable minerals 
indicates relatively rapid transport with little reworking (Maycock, 
1967), as does the jee of sedimentary, igneous (volcanic) and meta- 
morphic rock fragments. A western source was suggested. 

Approximately one hundred and sixty-feet below the top of the 
“Mannville in the southern part of the study area, a major coal seam 
occurs. This seam was correlated with seam 8 in the northern part 


of the area on the basis of its stratigraphic position, using the top 
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of the Mannville as the datum (figure 4.2). In both areas, the seam 
shows similar geophysical Jog characteristics with a relatively high 
resistivity and a negative SP deflection. Therefore, coal seam 8 at 
the top of the Borradaile Member in the northern area has its correla- 
tive in the south, although it is absent between Townships 36 and 39 


(figure 4.6). 


In the southern portion of the study area where Mississippian rocks 
crop out below the Mannville, the Deville Member is the basal unit of 
the Mannviile Group (figure 4.8). 

The Pense Formation, together with the underlying Cantuar Formation, 
forms the "basal Cretaceous sandy group" of Price (1963, p. 26). Accord- 
ing to Price and Ball (1971), the Pense Formation is present over most 
of the northern Williston Basin and has an average thickness of about 
150 feet. In the Duval potash shaft (Twp. 36, Rge. 6 W3M), the forma- 
tion is less than 50 feet thick (Price and Ball, 1971) and thins to zero 
a short distance to the northwest towards the area of Maycock's (1967) 
Undifferentiated Mannville Group. In the Duval shaft, the Pense Forma- 
tion is composed mainly of quartzose sandstone with numerous shale in- 


terbeds in the basal third of the interval. 


D. SOUTHWESTERN SASKATCHEWAN AND THE STUDY AREA 
In southwestern Saskatchewan, Christopher (1974) divided the Mann- 
ville Group into the Success, Cantuar and Pense Formations (figure 4.1). 
The Success Formation is composed of kaolin-indurated quartz sandstone 
with accessory chert. The Success Formation, plus the basal McCloud 


Member of the Cantuar Formation, are equivalent to the Lower Mannville 


of Alberta. 
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The McCloud Member of the Cantuar Formation is composed "... 
largely of autochthonous sandstones overlain by an estuarine to marine 
shale that reaches southward from east-central Alberta to the Inter- 
national Border in Saskatchewan as ria-like fingers of the marine [sic] 
Ellerslie Formation." (Christopher, 1975, p. 523). It is correlated 
with the Dina Member and the Cummings shale unit (Christopher, 1974, p. 49). 

In the present study area, deposition during McMurray time (Dina 
Member) as well as during Clearwater time (Cummings Member) was con- 
fined to north of Township 40 (figures 4.3 and 4.4). Based on the data 
available, the writer does not agree with Christopher (1975, text-fig- 
ure 4) that the Dina Member and the Cummings Member or their equivalents 
extend further south than Township 40. In this figure, Christopher 
(1975) shows the McCloud Member, in "valley form", extending southward 
from the "marine Ellerslie Formation" (Christopher, 1975, p. 523). The 
present writer could not recognize either the Dina or the Cummings (or 
their equivalents) on the geophysical logs in areas other than that 
shown in figures 4.3 and 4.4. Moreover, the "Ellerslie Formation" js 
probably not of marine origin (Williams, 1963, p. 356). 

Maycock (1967, p. 52), on the basis of petrographic similarity be- 
tween his Continental Facies in southwestern Saskatchewan and the Grand 

Rapids Formation, concluded that in western Saskatchewan "... the lowest 
beds in the 'Undifferentiated Mannville Group’ are ... most probably 
Aptian and/or Lower Albian in age and the Continental Facies and Pense 
Formation Middle Albian". On the other hand, Christopher (1974, p. 50) 
’ stated that "The sandstones of the McCloud Member ... in the Tidewater 
Southey Crown [Lsd 4-29-22-28 W2M] well is [sic] placed in the Barre- 


mian, and the superjacent shale (calcareous member) is dated upper 
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Neocomian, i.e., presumably Upper Barremian. The Dimmock Creek (glauco- 
nitic member) is given an Upper Barremian to Albian (?) but possibly 
Aptian age." However, both writers realized the difficulties in assign- 
ing a definite age to the Lower Mannville sediments due to the lack of 
faunal and/or floral evidence. Unless further work is done in this area, 
the evidence presented leads the writer to conclude that during the depo- 
sition of the McCloud Member and possibly the Dimmock Creek Member in 
Christopher's study area, deposition of the Undifferentiated Mannville 
and/or the Continental Facies had not even begun. Therefore, equating 
the Ellerslie sandstone with the McCloud Member (in "valley form") as 
shown in Christopher's (1975, text-figure 4) figure appears to be ques- 
tionable (see also "Interpretation of the Vertical Variability Maps in 


Chapter 7). 
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Chapter 5 
DEPOSITIONAL ENVIRONMENT OF THE MANNVILLE GROUP 


A. INTRODUCTION 

A complete regional synopsis of the depositional environments of 
Cretaceous time has been given by Williams and Stelck (1975) and Jelet- 
zky (1971). 

Following a long erosional period during pre-Mannville time, the 
topography in the study area had a maximum apparent relief of approxi- 
mately four hundred feet when Mannville deposition began (figure 5.1), 
being high in the south where the resistant Mississippian rocks cropped 
out (figure 5.2) and more gentle in the north where Devonian beds were 
exposed. This apparent relief may be partly due to contemporaneous 
salt solution of the underlying Devonian evaporite, leading to collapse 
of the overlying Mannville beds and hence, increases sediment accumula-~ 
tion locally. The term ‘apparent relief' is used because salt collapses 
occurred both syn- and post-deposition of the Mannville sediments. On 
the basis of the difference in relief, and with the Battleford Arch 
(approximately Twp. 40) as a dividing line (see discussion below), the 
area can be divided into northern and southern portions for discussion 


purposes. 


B. NORTHERN AREA 
The Dina Member was deposited during early Cretaceous (Aptian to 
Early Albian) time. The source of the member’was suggested by Nauss 


(1945, p. 1614) as possibly being from the underlying Devonian lime- 
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stone. However, in central Alberta, the corresponding McMurray Forma- 
tion has a Precambrian source (Williams et al., 1962). The Dina Member 
filled in the topographic lows, reducing the topographic relief before 
the transgression of the Clearwater Sea in which the Cummings Member 
was deposited. The Battleford Arch set the southern limit for the 
Clearwater Sea and therefore the Cummings Member. Within the study 
area, both Dina and Cummings are sand and together they subdued the 
relief further. Seam 3 was deposited near the top of the Cummings in- 
dicating a temporary emergence at that time. Westward toward the 
trough in central Alberta, the Clearwater Sea transgressed approximately 
to the present vicinity of Calgary (Williams and Stelck, 1975, text- 
figure 2) and stayed much longer than in the northeast. 

No attempt has been made to subdivide the rest of the Mannville 
Group into individual members. It suffices to say here that north of 
the Battleford Arch, post-Cummings deposition was mainly continental, 
as was the case in the Vermilion area to the west (Nauss, 1945). The 
lack of seams 4, 5, 6 and 7 in the northeast (figures 5.3a and b) was 
probably due to non-deposition. Seam 8 occurs at the top of the Borra- 
daile Member ("Sparky" sand jn oil industry terminology). According to 
Nauss (1945), the Borradaile is a near shore, marine-reworked sand and 
the overlying O'Sullivan Member is a deltaic sand in the area to the 
west. Therefore, seam 8 represents the emergence of the Borradaile 
Member with development of widespread swampy conditions (see figure 
4.6). No younger seams were found in the area (seams 9 and 10 are of 
very local extent in the south). The absence of coal in the O'Sullivan 


Member is probably due to the increase in sediment supply, leading to 
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non-swampy dry continental environments. At the end of O'Sullivan time, 
however, the rate of subsidence must have increased, leading to the in- 
vasion of the late Early Cretaceous sea and the deposition of the Joli 


Fou Shale which ended the history of the Mannville Group. 


C. SOUTHERN AREA 

South of the Battleford Arch, depositional environments during 
Mannville time were different from those farther north. In this area, 
the entire Mannville Group is of continental origin and contains ten 
coal seams each of limited areal extent scattered vertically throughout 
it (figure 5.3). Price et al. (1971, p. 20) described the distribution 
of coal in the Cantuar Formation as follows: | 

"The thickest coal bed occurs in strata considered to be 

river channel and intertidal sediments and not in typical 

‘coal measures’. Channel filling and the rounding of some 

wood fragments suggest that this ‘coal’ has accumulated 

from transported wood such as that found in log jams along 

river banks or collections of driftwood in shallow bays and 

not from the classically conceived mat of swampy vegetation." 

This drift origin may explain the erratic distribution of the coai seams 
in the study area which is reflected in figure 5.3a. 

In describing the Cantuar Formation in east-central Saskatchewan, 
Price et al. (1971, p. 19) noted that "... the formation is distin- 
guishable eastward as far as the Manitoba Escarpment and weli into the 
southern part of the Williston Basin. Northwestward, towards the re- 
gion of North Battleford [Township 43], the beds are included in the 


Mannville Group", thus suggesting a facies change in the vicinity of 


| Township 43. 
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Lower Cretaceous sedimentation in the south began with the deposi- 
tion of the Deville Member. According to Christopher (1974, p. 104), 
the base level of the Deville basins were higher than the levels of the 
master streams and hence, the basins formed an internal drainage system. 
It was this centripetal drainage system that eroded the Mississippian 
rocks and produced the cherty Deville Member. South of the Deville 
basins, the Swift Current Platform sloped gently to the south where the 
Success Formation was deposited. Deposition of the Cantuar Formation 
proper was brought about by uplift of the northern flank of the Swift 
Current Platform in the vicinity of Township 23 (Christopher, 1974, 
p. 104, figure 45), bringing the base level of the master streams higher 
than that of the Deville basins. Consequently, the drainage system 
turned to the north, breaching the Deville basins and the Mississippian 
outcrops at several places. On the Devonian strata to the nortn, the 
streams joined a master stream of subsequent valleys opening west-north- 
westerly into Alberta, where they join similar valleys, such as the 
'St. Paul Channel' mapped by Williams (1963). Within the study area, 
subsidence and deposition were probably in equilibrium, resulting in 
a low lying alluvial plain where streams flowed slowly northwestward, 
emptying into the east-central Alberta Plains (see Nauss, 1945, p. 
1627) and depositing their load as deltas and related features. On 
the alluvial plain, coal was deposited. 

The major seam in the southern portion of the study area is seam 
8, which, however, is absent from Townships 37 to 39. In view of the 
uniform thickness of the Mannville sediments overlying seam 8 (figure 


4.2) and the fact that no unconformity was detected at the base of 
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these sediments, the writer speculates that the absence of seam 8 was 
due to non-deposition. 

Deposition of the Cantuar Formation in the south (approximately 
south of Twp. 32) was brought to a close by the shallow marine trans- 
gression which produced the sands of the Pense Formation. Within the 
present study area, the Undifferentiated Mannville Group (equivalent 
to the Cantuar Formation) is overlain by the late Lower Cretaceous Joli 
Fou Formation. 

In general, central Saskatchewan, during Lower Cretaceous time, 
was high. In the northern portion of the study area (north of Twp. 40), 
near shore and deltaic environments alternated through Mannville time. 
In western Saskatchewan (from approximately Rge. 24 to 28 W3M) along 
the provincial border, a fluvial plain environment dominated the area 
(Maycock, 1967), where coal seams were deposited (figures 5.3a and b). 
Across the border to the west the McMurray and the Grand Rapids Forma- 
tions were deposited in fluvial and deltaic environments while the me- 
dian Clearwater Formation is a marine shale. The eastern part of the 
map area remained above sea level most of the time with occasional ma- 
rine transgressive tongues and sparse coal deposition. 

After the deposition of seam 8, the resulting topography was con- 
siderably subdued, with an apparent relief of about 110 feet (figure 
5.4). The absence of seam 8 in Townships 37-39 is due to non-deposi- 
tion. Uniform deposition continued thereafter until the end of Mann- 
ville time. The top of the Mannville shows an apparent relief of 180 


feet (figure 5.5) as compared with 425 feet in the sub-Cretaceous ero- 


sional surface (figure 5.1). 
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D. SPECIAL FEATURES 

From the descriptions of the stratigraphy and the depositional en- 
vironments, three special features stand out clearly within the study 
area. They are as follows: 

(1) the southern limit of both the Dina and the Cummings Members 
is in the vicinity of Township 40 and trends approximately east-west, 
whereas the equivalent McMurray and Clearwater Formations in central 
Alberta extend much farther south; 

(2) the absence of seam 8 between Townships 37 and 39; 

(3) the completely different depositional environments between the 


northern (north of Twp. 40) and the southern parts of the area. 


E. THE BATTLEFORD ARCH AND THE MANNVILLE GROUP 

The Battleford Arch in central Saskatchewan is the northeastern 
extention of the Sweetgrass Arch in southern Alberta. In his discus- 
sion of basement contro! of Cretaceous sand sequences in western Canada, 
Stelck (1975) has shown the distribution of Precambrian basement fea- 
tures (figure 5.6). During earliest Cretaceous (Neocomian) time, ma- 
rine sediments (Monteith, Beattie Peaks and Monach Formations) were 
confined to the Keg River Low north of the Peace River Arch. South of 
the arch, deposition of the Blairmore Group did not start until Aptian 
time, and in the West Alberta Basin, the Blairmore lies unconformably 
on Jurassic or older strata. The basal conglomerate (Cadomin Formation) 
of the Blairmore Formation only exists in the Foothills area and jis 
~ absent under the Plains. The Peace River Arch apparently was an area 


of erosion (see Stelck, 1975, p. 432, text-figure 4). In mid-Albian 
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time, the arch became a demarcation line for facies change; continental 
and near shore deposits in the south gave way to marine facies in the 
north. 

A similar situation occurred in the present study area. Deposi-~ 
tion of the Dina Member began during Aptian time and its distribution 
was confined to the area north of Township 40. Both in the study area 
and in Nauss' (1945) area, the member has a sandy lithology and thick- 
ens towards the north. After deposition of the Dina Member, the Clear- 
water Sea deposited the Cummings shale in east-central Alberta and sand 
fs the study area. Like the Dina Member, the Cummings does not extend 
farther south than Township 40. Therefore, deposition during Aptian to 
Early Albian time was confined to the north of the study area. The 
barrier that limited the southern extension of the two members is at- 
tributed by the writer to the presence of the Battleford Arch. 

Deposition of the continental Cantuar Formation in the southern 
part of the study area did not start until late Clearwater time when 
the sea regressed to the north. The overlying Islay, Tovell and Bor- 
radaile Members are beach, deltaic and reworked beach sand deposits 
respectively, whereas the Cantuar Formation remained continental. 

This difference in Pcie: change, though not as striking as in the 
Peace River Arch example, is significant enough to indicate the pres- 
ence of a demarcation barrier. Coal seam 8 ("Wainwright coal", Wilson, 
1947) at the top of the Borradaile Member is absent between Townships 
37 and 39 due to non-deposition, a fact which also suggests a high in 


this vicinity which can be attributed to the presence of the Battle- 


ford Arch (figure 4.2). 
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As stated by Stelck (1975), basement control over the deposition 
of Lower Cretaceous sand is subtle. The Battleford Arch was high dur- 
ing the deposition of the Mannville Group in the study area and con- 


trolled the distribution of the Dina and Cummings Members and seam 8. 
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Chapter 6 
STRUCTURAL GEOLOGY 


The structure of the Mannville Group, summarized on figures 6.1, 
6.2 and 6.4, essentially reflects the topography of the underlying 
sub-Cretaceous erosional surface. 

Figure 6.1 shows the present configuration of the sub-Cretaceous 
surface. The major dominant feature is the sinusoidal positive feature, 
named here as the Unity Ridge. Other ridges and valleys run subparailel 
to this main ridge. These topographic high and low features of the land 
surface are the results of erosional effect. 

Remnants of the Mississippian rocks create east-west trending ridges 
in the south (Maycock, 1967, p. 55). The south-trending valley in the 
south of the map area is probably one of the northern arms of the Plato 
Valley named by Christopher (1975, text-figure 4). 

Structure on the top of the Mannville (figure 6.2) is different from 
that of the base. The sinusoidal Unity Ridge may be seen in the north 
only, reflected in two offset positive features. Since no major fault 
was found in the area, the offset is probably a depositional feature and 
not tectonic. An indistinct, north-trending nose is present in the 
southeast corner of the map. In general, the features which are very 
distinct in figure 6.1, are covered up by the Mannville sediments and 
only the barest outlines remain. The surface has acquired a much more 
uniform southwesterly dip. 

Depressions E and F in the northeast corner of the area show per- 


sistent closure of up to 250 feet through the entire Lower Cretaceous 
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(figures 6.1, 6.2 and 6.3). They developed in post-Early Cretaceous time 
Since isopach maps of the Mannville Group (figure 6.4) and from the Base 
of Fish Scales to the top of the Mannville (figure 5.5) do not show any 
thickening in that area. The depressions are attributed to the collapse 
of overlying strata due to the removal of underlying Devonian salt. 

Depression G in the south of the map area decreased in areal extent 
through Mannville time and it persisted to at least the end of Early 
Cretaceous (figure 6.3). Its relief at the base of the Mannvilie (figure 
6.1) is approximately 50 feet. The isopach map of the Mannville Group 
(figure 6.4) shows only a 25-foot thickening in the same general area, 
indicating salt collapse started during Mannville time and continued through 
the Early Cretaceous. Similarly, depression H has a small areal extent 
with a relief of less than 50 feet at the base of the Mannville (figure 
6.1). The isopach map of the Mannville Group (figure 6.4) also shows less 
than 50 feet of thickening of the sediments in the same area. Syn-depo- 
sitional salt collapse was the origin of this depression which was fiiled 
in during Mannville time so that it does not show up in the structural map 
of the top of the Mannville (figure 6.2). 

The structural section in figure 6.5 shows the regional southerly 
tilt of the map area, compared with the stratigraphic section (figure 4.2). 
In discussing the intermediate structures in south-central Alberta, 
Robinson et al. (1969, p. 2367) noted that "... the intermediate structures 
of the Cretaceous surfaces probably resulted not from differential compac- 
tion but from tectonic movements that continued throughout the Cretaceous 


| and possibly into Early Tertiary". This may also have been the case in 


west-central Saskatchewan. 
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Chapter 7 
COAL RESOURCE ESTIMATION 


A. FACTORS AFFECTING RESOURCE ESTIMATION 

Several factors have a direct effect on the estimate of the coal 
resources of the study area. These factors are summarized below and 
the reader should keep them in mind when considering the tonnage of coal 
arrived at in the estimation process. 

Of the 242 townships in the area, 31 contain no logged boreholes. 
These townships occur mostly in the central part of the study area. 
Since resource estimation is calculated on a one-well-per-township basis, 
these 31 townships are automatically exciuded from the calculation - 
that is, contribute zero value. The total coal thickness in each town- 
ship, plus the areas with no data, are shown in figure 7.1. From the 
formula for resource estimation (see below), it may be seen that each 
addition of a one-foot seam per township will increase the estimate by an 
increment of 40.7 million tons, hence, the lack of data for the 31 town- 
ships may well affect the estimate significantly. On the other hand, 
because of the laterally discontinuous character of the seams (see sec- 
tion D, this chapter), the extrapolation of data from a single well into 
a zone of influence of a township involves a large degree of uncertainty 
and therefore more control points would undoubtedly improve the result. 

In 53 townships in the northern part of the study area, the deepest 
available boreholes did not penetrate the entire Mannville Group. Most 
of these holes were drilled only to the Tovell Member ("Rex sand" in oi] 
industry terminology). As mentioned in the discussion on the strati- 


graphy of the northern part of the area, a major coal seam (seam 3) occurs 
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Fig 7.1 TOTAL COAL THICKNESS DISTRIBUTION. 
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near the top of Cummings Member in the north. Hence, the final estimate 
is subject to revision should additional data from deeper boreholes be- 
come available. 

The frequency of seam occurrences is shown in figure 7.2. The two 
major seams, 3 and 8, were encountered 22 and 59 times, respectively. 
Seam 6, which underlies the "North Hoosier Basal Blairmore Sand Pool" 
(White, 1974, p. 128), was found 17 times in the 236 wells. 

The geophysical logs used in this study have a resolution of approx- 
imately 2 feet and seam thicknesses were read to the nearest foot. How- 
ever, because of the low resolution of the logs, seams less than 2 feet 
thick were not identified. In the study area seam thicknesses range from 
2 to 9 feet. The mean value is 3.9 feet with a standard deviation of 1.8 
feet. Most seams are either 2 feet or 4 feet in thickness (figure 7.3). 

- As discussed in Chapter 3, each coal seam as determined from logs 
was assigned a reliability index of 1, 3 or 5 (also see Yurko, 1976, 
pp. 11-12). Figure 7.4 shows that 78.1% of the seam identifications were 
based a a resistivity log and a sonic log (reliability index 1), 18.8% 
on a resistivity log only (reliant ites index 5), and 3.1% on resistivity 


and neutron logs (reliability index 3). 


Figure 7.5 shows the total coal thickness map plotted from one 
hundred and twenty-eight data points including reliability indexes 
1, 3 and 5. If only reliability index 1 is used, only one hundred 
data points are included, but the isopach map shows a very similar 
pattern. Both maps show a general trend of thickening toward the 


west of the study area. 
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The isopach map of seam 8 is shown in figure 7.6. This map includes 
fifty-nine data points with reliability indexes 1 (41 data points), 3 
(2 data points) and 5 (16 data points). Again, when only reliability 
| indexes 1 and 3 are used, the isopach map shows a similar trend of 


thickening towards the Saskatcnewan-Alberta border. 


B. METHODS OF ESTIMATION AND RESULTS 


The estimate of the coal resources in the study area is based on the 


following formula (Averitt, 1969): 


i SOaX eV OL EXD xe G 
Tons per township = TTD COOe a 


where 36 is number of sections per township; VOL is the volume of coal 
in cubic feet in one section (one square mile), ji.e., 5,280 x 5,280 x 
thickness in feet; DW is density of water (62.4 1b/cu ft); and SG is 
specific gravity of coal (1.30 gm, an average selected for the subbit- 
uminous coal in Saskatchewan). 

Substituting in the formula, one obtains a value of 40.7 million 
short tons per township for a one-foot coal seam. 

An estimated value of 19,400 million short tons of coal (including 
all seams 2 feet thick or thicker) was obtained for the Mannville Group 
in the study area. The major concentration of coal is near the Alberta- 
Saskatchewan border from Townships 31 to 35 and 39 to 43 (figure 7.1). 
If only seams of 3 feet thick or thicker are considered, the estimated 
value drops to 16,800 million tons, i.e., 86.6% of the initial value. 

If only seams of 5 feet or greater thickness are considered, the estimated 


value is further reduced to 8,500 million tons or 43.8% of the total. 
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Table 7.1 summarizes these values. The resource estimates on the two 
major coal seams (seams 3 and 8) are shown in Tables 7.2 and 7.3. If 
Seams of 2 feet thick or thicker are considered, seam 3 constitutes 16.1% 
of the total resource estimated for the study area while seam 8 consti- 
tutes 49.3%. 

Latour and Chrismas (1970) estimated a total of 12,014 million tons 
of coal resources in place in southern Saskatchewan. This coal is of 
Tertiary and/or Late Cretaceous age and includes all seams greater than 
three feet thick. Within the study area, the corresponding resources in 
the Mannville Group are 13,879 million tons. 


Yurko (1976) made a similar estimate of the coal resources of the 


Cretaceous System in Alberta based on one petroleum borehole per township. 


Table 7.1 


Resource Estimate of Mannville Coal 
(in millions of short tons) 


Seam 
Thickness 
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Table 7.2 


Resource Estimate of Seam 3 
(in millions of short tons) 


Seam Percent 
Thickness of Total 


Jabless73 


Resource Estimate of Seam 8 
(in millions of short tons) 


Percent 
of Total 


Seam 
Thickness 


100.0% 
85.6% 
74.2% 
47 .0% 
34.3% 


Table 7.4 


Coal Resources in Various Areas 
(in millions of short tons) 


Mannville Group 


(Lower Cretaceous) Tertiary and/or 


Late Cretaceous 
Alberta | Southern Sask. 
Seam Yurko Latour and 


Thickness (1976) Chrismas (1970) 


627,978 - 
12,014 


401,724 - 
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He obtained a total estimate for the entire Mannville Group over the whole 
province of 627,978 million tons for seams 2 feet thick or thicker; his 
estimate for seams thicker than 5 feet was 401,724 million tons which is 


64% of the total (table 7.4). 


C. VERTICAL VARIABILITY OF THE COAL SEAMS 

In the Mannville Group, up to 10 seams are distributed through 
up to 280 feet of strata. A study of the vertical distribution of coal 
seams in such a stratigraphic section can provide valuable insights into 
palaeoenvironmental conditions in Mannville time. The technique employed 
involved calculation of the centre of gravity of coal seams in each 
borehole using the procedure of Krumbein and Libby (1957), and presenting 
the data in map form (figure 7.7). 

The centre of gravity map shows the weighted mean position of the 
coal in terms of its distance from the top of the Mannville Group. I[t 
should be noted that the data have been extrapolated into areas where no 
coal or no control points occur. Within the study area, the centre of 
gravity ranges from 140 to 380 feet below the top of the Mannville. To 
express the centre of gravity in terms of the percentage of the thick- 
ness of Mannville as measured from the top of the group, the relative 
centre of gravity is plotted (figure 7.8). In the northern part of the 
area, the relative centre of gravity falls mainly in the 40% to 60% 
range, indicating that most of the coal is approximately at the midpoint 
of the Mannville Group. In the southern area, the relative centre of 
gravity drops to the 60% to 80% range, showing that the seams are aggre- 


gated near the base of the Mannville (figuinea. 2). 
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Information on how the seams are actually distributed in the Mann- 
ville Group is provided by the standard deviation maps (figures 7.10 and 
7.11). The standard deviation, in this case, is that thickness of the 
formation, measured from the centre of gravity, which includes 68% of 
the coal. The relative standard deviation map (figure 7.11) shows the 
percentage of the Mannville strata being included in one standard devia- 
tion. Zero standard deviation occurs where there are no control points, 
or only one seam is present in the well. As the value increases, the 
Spreading out of the seams also increases. The major portion of the 
coal-bearing area lies below the 60-foot contour, indicating that within 
this area, 63% of the coal occurs within 60 feet of the centre of gravity 
(figure 7.10). However, in the southwestern part of the area, plus or 
minus 60 feet (i.e., 120 feet in total) is almost half the thickness of 
the entire Mannville Group. This can be illustrated by reference to the 
relative standard deviation map (figure 7.11) which shows the percentage 
of the Mannville Group included in one standard deviation. Like standard 
deviation, zero relative standard deviation indicates no coal or only one 
seam in a well. As the value increases, the percentage of the formation 
included in a standard deviation also increases. For example, within the 
20% contour, in order to account for 68% of the coal in a well, 20% of 
the formation on either side of tne centre of gravity must be included 
that is, a total of 40% of the formation. Two standard deviations will 
account for 95% of the coal, therefore inasmuch as the major part of the 
study area lies below the 20% relative standard deviation contour, on the 
~ average, 95% of the coal in a well occurs within an 80% thickness inter- 
val of the Mannville sediments. It is apparent therefore, that the coal 
seams in the Mannville Group are scattered throughout the group and are 


not aggregated in any one part of it. 
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D. INTERPRETATION OF THE VERTICAL VARIABILITY MAPS 


A comparison among the total coal thickness distribution (figure 
7.1), the centre of gravity map (figure 7.7), the standard deviation 
map (figure 7.10) and the isopach map of the Mannville Group (figure 6.4) 
shows that where major coal accumulation occurs, the following conditions 
generally exist: 

(1) The Mannville sediments are thick, indicating a palaeotopographic 
low (figure 6.4); 

(2) The number of seams and/or the net thickness of coal (figure 7.1) 
has a direct relationship to the centre of gravity of coal (figure 7.7). 
When either the number of seams or the net coal thickness is large, the 
value for the centre of gravity is also large, ji.e., the centre of gravity 
is further down from the top of the Mannville; by the same token, when 

(3) The number of seams and/or the net thickness of coal (figure 7.1) 
is large, the standard deviation of coal is large, i.e., the vertical 
scattering of the seams increases. 

It has to be emphasized that these conditions are generalized rela- 
tionships. Statistical treatment af the geophysical log data using linear 
regression shows that ein is little correlation among the parameters 
discussed. However, by superimposing the maps, the qualitative relation- 
ships mentioned can be shown. 

These conditions can be explained by hypothesizing the depositional 
environments during Mannville time in the study area as follows. The 
palaeotopographic low areas became the sites of either swampy areas or 
| stream channels. In these swampy areas, coal was deposited. Log jams 
along river banks (see Price, 1971) may also have been a source of coal 


accumulation. As sedimentation of the Mannville Group continued, the 
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loci of these low-lying areas migrated laterally, leading to the present 
areal distribution pattern of the coal. 

Condition 2 arises since the majority of the seams (seam 1 to 8) 
were deposited in the lower portion of the Mannville Group. This indi- 
cates that the supply of sediments was in pace with subsidence, leading 
to equilibrium (or stable) conditions for the continual deposition of 
coal. However, the decrease in coal deposition after seam 8 signifies 
either an increase in sediments supply or a rapid subsidence of the ba- 
sin. The former case gives rise to a dry continental environment, where- 
as rapid subsidence will lead to marine transgression. Since the portion 
of the Mannville Group above seam 8 is of continental origin, tne writer 
adopts the theory that sediment supply increased after the deposition of 
seam 8. 

If the loci of deposition of coal was static for a long period of 
time during the deposition of the Mannville Group, there would be fewer 
seams and each seam would probably be much thicker. However, this was not 
the case within the study area as indicated by the standard deviation map 
(figure 7.10). and the net coal thickness map (figure 7.1). The frequent 
change of depositional loci led to increasing scatter, both geographically 
and stratigraphically, of coal distribution in the Mannville Group. 

A study of the position of the centre of gravity and the standard 
deviation of coal (see also figure 7.9) in the study area further shows 
that north of the Battleford Arch, the Mannville sediments can be divided 
into three zones; a coal-bearing zone being sandwiched between lower and 
upper non-coal-bearing zones. The lower non-coal-bearing zone consists 
of the Dina and the Cummings Members, the overlying coal-bearing zone 


extends approximately from the top of the Cummings Member to the top of 
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the Borradaile Member (i.e., seam 8), and the upper non-coal-bearing zone 
includes those Mannville beds overlying seam 8. 

South of the Battleford Arch, the lower non-coal-bearing zone is 
absent and the median coal-bearing zone is bounded by the sub-Cretaceous 
unconformity below and seam 8 above. The overlying non-coal-bearing zone 
has the same vertical extent as its northern counterpart (figure 7.9). 
The correlation between the coal-bearing zones both north and south of 
the Battleford Arch further supports the hypothesis that during the depo- 
sition of the Dina and the Cummings Members in the northern part of the 
study area, deposition of the Undifferentiated Mannville Group in the 


south had not begun (see discussion in section D, Chapter 4). 
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Chanter 8 
AREAL VARIATION OF THICKNESS OF SEAM 8 


A. THEORY AND APPLICATION 

Matheron (1963) defined a regionalized variable as one that has a 
definite value at each point in space. Examples for this kind of vari- 
able are coal seam thickness, ore grade, water chemistry in an aquifer, 
etc.. The value of these variables may change from one point to another 
and the spatial relationship of this change (i.e., how rapid is the 
change from point to point) is important to the geologist. The tool to 
examine this spatial relationship is the variogram which is a statisti- 
cal measure of the degree of areal variation or lateral continuity of 
a variable. A full mathematical treatment of the variogram is outside 
the scope of this thesis. The reader is referred to the publications by 
Matheron (1963), David (1974) and Huijbregts (1975). The following is 
an account of the basic concept of the variogram. 

Suppose one wants to compare, say, coal seam thickness at two points 
(x) and (xth) at a distance h units apart, one would examine the thick- 
ness difference at the two points without being concerned with the sign 
(positive or negative) of the difference. The simplest way to eliminate 
the sign is to square the difference. To express this in mathematical 


form: 
2y(h) = [th(x) - th(xth)] 


where 2y(h) = the variogram; th(x) = thickness at point x; and th(xth) 
= thickness at point (xth), h units away from point x. Usually, the 


— 


semi-variogram, y(h), is used. 
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The variogram, 2y(h), is a vector and takes the direction of h, the 
separation of the points, into account. Variogram values refer to pairs 
of points; for example, if one has four data points A, B, C, and D, one 
will have a total of 6 variogram values, for the pairs AB, AC, AD, BC, 
BD, CD, that is, [(n-1)+(n-2)+...+(n-n)] values. To simplify matters, 
variogram values are grouped into classes according to the distance be- 
tween the pairs and the average value of the class is taken as the repre- 
sentative value for the class. Hence, for a particular class, say within 


0 to 500 feet, the variogram value can be written as: 
[th(x) - th(xth)]? 


where (h)-denotes the direction of the variogram value, N is the total 
number of data point pairs in that class, j.e., 


2(h) = average [th(x) - th(xth) ] 


When semi-variogram values are plotted versus distance between points 
four general aces OG curves (Matheron, 1963) may be recognized (figure 
8.1). The continuity of the variable under consideration is reflected by 
the rate of increase of y(h) for small increments in distance. Type A has 
a parabolic trend at the origin and the variable under consideration (e.g., 
bed thickness) has a high continuity. Type B shows that the change of y(h) 
is uniform, typifying, for example, a relatively uniform sedimentary de- 
posit where changes usually occur very slowly. In other cases, again using 


a sedimentary example, deposition may be in the form of lenses or even 
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nodules, often concentrated in small areas, with rapid changes occurring 
over very short distances. Curve C shows the effect of this rapid change; 
the positive intercept (Co) is known as the nugget effect. The nugget 
effect may occur for a number of reasons, such as poor sampling resolution 
(for example of geophysical borehole logs) or erratic distribution of the 
variables. It increases the variance of the variable under consideration 
by the amount equal to Co and also implies that if one goes back to-sample 
the same site, the results will differ (see David, 1974, p. 62). When the 
semi-variogram values, y(h), are independent of each other regardless of 
the distance between two samples, a curve of type D is produced. This is 
common in gold deposits for example and, it is called the random (or pure 
nugget effect) new 

The zone (or range) of influence of a data point is the distance be- 
yond which the influence of the data point disappears. In the variogram, 
this zone is represented by the distance at which the curve reaches a 
plateau (the sill). Along the sill, the rate of change of the variogram 
is independent of the distance. In curve D, the zone of influence is zero 
for all practical purposes. 

As has been mentioned, the variogram, 2(h), is a vector, and there- 
fore when specifying the zone (or range) of influence, the direction in 
which the zone is effective must be taken into consideration. For this 
reason, the word "zone" (or range) is preferable to "radius" which implies 
the influence is the same in all directions. Because of this directional 
component, if one compares variograms taken in different directions and 
plots the corresponding results, the zone of influence for the different 


curves may not be the same. If the zone of influence is the same in all 
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directions, the condition is called isotropic, otherwise the term anis- 
tropic is applied to the condition where the zone of influence varies 
with direction. In figure 8.2, for example, the zone (or range) of in- 
fluence in the north-south direction is approximately twice as much as 
that in the east-west direction. By using the zone of influence and its 
direction, one can reconstruct the major axis of continuity of the de- 
posit with some degree of certainty. Venter (1975) has discussed the 
use of the variogram to increase the reliability and validity of the 
conventional polygon method of calculating coal reserves. This is 
achieved first by limiting the maximum area-of-influence by any one 
drill hole, and then by assigning the width of the boundary zone of a 
deposit according to the zone of influence obtained from the variogram. 

In the present study, it is intended to use the variogram to gain 
some insight on the palaeodrainage direction during Mannville time in 
the study area. It is hypothesized that if coal was deposited on the 
stream floodplains, deposits would have a greater extent and consistency 
parallel to the stream flow than in the transverse direction. If such 
is the case, then the continuity of the coal seam thickness will be 
greater in the direction parallel to the stream flow. This continuity 
should be reflected by the zone of influence in the variogram; the greater 
the continuity, the longer the zone of influence. Hence, by obtaining 
the zone of influence in various directions, it may be possible to recon- 
struct the palaeodrainage direction. 

One hundred and one wells (figure 8.3) were used to depict the 
thickness of seam 8 in part of the study area. The resulting variograms 


are presented in figures 8.4 to 8.10. Seam 8 was chosen for the vario- 
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VARIOGRAM WITH 2B PIELD OF 90. DEGREES IN EACH DIRECTION 


STEP IN METERS @ 500.0000 
UPPEB LIAIT PCB Z ® 30.0000 
CEMERAL AZAW OF 2 ; 8 3.9208 
GENERAL WARIANCE OP 2 8 0.9937 
DISTANCE IN METERS UB. OP PAIRS DRIFT VARIOGRAN AYEBAGE DISTANCE 
0 ---= 500 72 - 0.097 27831 372.3 
500 ---- 1090 172 0.186 0.7616 775.3 
1CCO ---= 1506 217 =0.005 0.6705 1266.2 
1500 ---= 2690 : 253 0.170 0.6739 1745.6 
2090 ---= 2590 264 0.987 0.9375 2253.7 
2590 ---= 3000 228 0.083 0.8487 2751.2 
3C0C0 -=-- 3570 177 0.102 0.9096 . 3299.2 
3500 ---- 4090 168 -0.131 0.8278 3730.3 
809) —--"8500 139 -0.098 0.9532 6237.1. 
8509 ---= 5000 195 9.021 6.8247 5755.2 
50900 ---- 5500 138 -0.355 1.1268 5278.0 
5590 ---- 6C00 186 0.110 1.0685 5742.0 
6000 =--= 6590 182 -0.220 1.2747 6238.6 
6596 ---= 7000 172 0.453 - 1.4593 6732.6 
7000 —-=-- 7500 179 0.011 1.0726 7238.1 
75€0 ---= 800¢ 123 -0.317 1.3618 7737.2 
8000 ---- 8500 80 -0.875 1.5590 8213.2 
8500 --== 9000 32 -0.375 4.1562 8698.6 
8000 =<--- 9500 9 -0.333 0.8333 9178,2 
— 9560 ----10000 ? 2.000 2.0000 9526.6 
ah) 10600 —---10590 1 0.0 0.0 1003527 
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1.430 © - 
1.399 @ 
1.369 @ A 
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1.278 ® 
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4.217 © 
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1.156 ® 
1.125 @ © 
1.035 ® A 
1.065 ® s e 
1.038 @ 
1.008 © 
a aaa aa EDM TST EST 7 STREETS ISS OEE EO ERS, 
0.9483 © ® variance 
0.913 ® © 
C.882 © bg 
0.852 * 
¢.821 ® * ~ = 
0.791 © 
0.760 © e 
0.730 © 8 
0.700 © 
0.669 © * - 
0.639 © 
0.608 ® 
0.578 * 
0.548 © 
0.517 © 
C.ue7 © 
9.459 ® 
0.426 © 
0.395 © 
0.365 @ 
0.335 © 
0.308 # 
0.278 ® 
C.243 © 
¢.213 ¢ 
0.183 © 
¢.152 © 
0.122 @ , 
0.091 © 
0.061 © 
0.630 © 
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STEP I= METERS 


OPPER LIMIT FOR Z 


GEWEBAL AZAM OP Z 


16899 
7.460 
1.430 
1.399 
1.369 
1.338 
1.308 
1.278 
1.287 
1.217 
1.186 
1.156 
42.125 
1.095 
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BITA A PIPLD oF 90. DEGREES IN EACH DIRECTION 


VARIOGRAM 

# 500.0000 

® .30.00C0 

§ 3.9208 

GENERAL VARIANCE OF Z  ] 0.9937 
DISTANCE IN METERS NB. OF PAIRS DRIFT 
C =--- 500 73 0.055 
$00 ---= 1€00 172 0.279 
10690 ---- 15C€¢ 228 0.CAC 
1506 =--= 260¢ 248 0.028 
2000 e--= 25099 297 0.030 
25C0 =--= 3090 285 0.172 
3600 -=== 3590 232 C.192 
3500 <== 2000 238 0.976 
BCOC <== 4500 210 0.019 
§5006 ---- 5000 192 0.193 
5020 === 5500 155 “0.071 
550C === 6€00 126 0.071 
6000 ---= 6500 133 90.258 
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$251.7 
5732.8 
6227.8 
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® 
VARICGRAM WITH A PYELD OF 80. DEGRRES IR FACH DIRECTIOV 


STEP IB RETERS @ 500.0000 
UPPER LIAIT FOR ¥ i @ 30.0000 
GEOFRAL BEAM OF Z ) 3.9208 
GEEPRAL VARIANCE OF 2 6 0.9937 
DISTANCE I® RETERS BB. OF PAIRS DRIPT VARIOGRAR AVERAGE DISTANCE 
0 -e-= 509 Ta ~0.081 0.9595 372. 
500 -=-= 1006 178 0.186 6.8708 782.9 
iCOO s--- 1500 209 0.038 0.8568 1268.8 
1500 -=-== 2000 218 -0.061 1.1238 1753.2 
2690 ---- 2595 279 =0.129 1.200? 2236.7 
25(0 ---= 3¢00 2 259 6.228 en tiow: 2731.1 
3000 s--= 3590 233 -0.008 1.0393 3238.9 
359C ---= 4060 293 0.072 0.858C 3736.6 
8690 e--= Ls0G 204 0.162 C.6867 4251.3 
8500 ---= 5000 786 0.164 0.6778 ‘ 8750.8 
5CC0 -2-=- 5569 135 0.108 0.9556 5242.2 
5500 --== 6000 106 0.396 1.1792 : 5730.9 
60900 ---- 6500 68 0.152 0.9779 6193.8 
6560 --== 7060 4&2 -0.500 6.9465 6799.5 
7000 -=-=- 7590 22 0.182 0.9545 7207.8 
500 ---= 8696 18 -0.278 0.8056 7719.14 
600G =-==- 8590 17 -0.298 0.9706 8230.2 
8500 ---= 9000 8 -0.750 1.75060 8758.38 
9000 ---= 9500 6 0.0 6.8333 9199.3 
$506 =--=-10000 1 -2.000 2.0606 9526.6 
¥ (h) 10000 --=~10500 4 0.0 0.0 10035.7 
4 1.899 # 
1.469 © 
1.430 © 
1.399 ® 
1.369 © 
1.338 © 
1.308 © 
1.278 ®@ 
1.287 8 
1.217 © 
1.186 ® ® 
1.1755 © iS 7 
41.125 ® 
7.695 ® fe 
7.065 © F f 
1.039 © Sd Variance 
1.008 ee ne ee 
0.973 ® $ = 
0.963 & © x 
C.913 ® 
0.342 & om 
0.852 * s s 
0.821 @ 4 
0.791 © 
0.760 ® 
0.730 @ 
C.703 &, : 
0.669 ® Ms ee 
0.635 ® 
0.608 ®@ 
0.578 ¢ 
0.548 ®& 
C.517 © 
C.987 © 
0.456 © * 
0.826 ® 
0.395 ® 
2265 @ 
C.335 ® 
0.308 © 
0.278 © 
0.283 ® 
0.213 © 
0.183 © 
0.152 * 
0.122 # 
0.091 © 
0.6613 © 
0.630 ©. 
0.000 ®* 
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GENERAL VARIANCE OP Z 


¥ (h) 
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1.299 
1.460 
1.339 
1.399 
1.369 
1.338 
41.308 
4.278 
1,287 
1.217 
1.1386 
1.158 
4.125 
1.095 
1.065 
1.038 
1.008 
0.973 
0.943 
0.913 
0.832 
0.852 
0.821 
0.791 
0.760 
0.730 
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VARIOGRAM 

§ 500.0000 

6 30.0009 

6 3.9208 

\] 0.9937 
BB. OP 
590 7S 
1¢00 176 
18¢90 202 
2c00 2Cco 
2500 219 
3000 207 
3500 213 
“too 215 
4sco 200 
$c00 164 
5590 128 
6C00 87 
65C0 36 
7000 21 
7500 7 
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s 3 
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PRIPT 


°0.213 
0.057 
-0.005 
-0.C€95 
-0.278 
C.188 
-0.047 
0.028 
-0.160 
0.122 
-0.055 
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0.000 
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YARICGRAS AVERAGE DISTANCE 
2. AKGT 367.5 
1.0304 766.9 
1.0965 1250.5 
1.0675 1757.0 
1.2968 2283.7 
1.1232 2735.9 
0.9906 3233.2 
0.8419 3780.9 
0.7409 8262.8 
0.8649 3740.3 
1.0430 5218.9 
1.4195 57C7.6 
1.1389 6195.7 
1.0716 6731.3 
1.2857 7181.9 
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gram study because the area has very good well denisty and the seam is 
relatively widespread, some of the wells have a complete suite of logs 
and most of the boreholes penetrate the seam. More data are available 


than for any other seam in any other area. 


B. RESULTS AND INTERPRETATION OF THE VARIOGRAMS 

In statistics, the variance of a variable (e.g., seam thickness) 
is a measure of how the variable is spread out around the mean. The 
more commonly used value is the standard deviation which is the square 
root of the variance. In this variogram study, three kinds of variance 
(figure 8.11) are involved and are discussed below. 

1. Variance of sampling error is represented by the nugget effect 
(Co) in the variogram curve (David, 1974, p. 62). In this study, all 
the variograms show a high nugget effect value meaning that in determin- 
jing the seam thickness from the geophysical logs used, an error was in- 
troduced by the sampler. This error arises because the resolution of 
the logs used is about two feet and the thickness was read to the nearest 
foot. With better log resolution, the accuracy could improve, hence 
decreasing the nugget effect. 

2. Variance due to natural random variation of the seam thickness 


is here called sample variance. This variance can be calculated by the 


formula: 


oe n(n x th)2 - (sn x th)? 
ig N 


where §2 = variance, n = number of observations that have thickness th, 


th = thickness of seam, and N = total number of observations. 
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In this study, the sample variance is 0.9937 feet®. However, this 
sample variance is directly affected by the sampling error. If an error 
is introduced during determination of seam thickness, the sample vari- 
ance calculated becomes an unreliable estimation of the natural random 
variation of the seam thickness. In view of the high nugget effect 
value (approximately 0.7 to 0.8 ft?) obtained in this study, 0.9937 feet2 
is probably not representative of the true vaiue of the sample variance. 
3. The third type of variance is that of the error made when esti- 
mating, say, seam thickness at point (xth) using the thickness at point 
X,; which is h units from the former point. This variance is known as 
estimation variance. Huijbregts (1975, p. 41) has shown that the vario- 


gram being the average quadratic difference between thickness at (xth) 


and <thiteknesswat x, Tie 23 


2v(h) = average [th(x) - th(x+h) 1 
this function is an estimation variance. The estimation variance in all 
types of variograms (figure 8.1), except the random type, increases as 
the distance between the pair of points increases. It means that in ex- 
trapolating a data value, the farther away from the contro! point, the 
higher is the variance introduced into the estimated value. The maximum 
of the estimation variance is represented by the sill of the curve. 

Along the sill, the estimation variance is independent of the distance 
between the point-pair. The zone of influence is defined as the distance 
at which the variogram curve reaches a platedu (the sill), i.e., where 


the estimation variance reaches the maximum. 
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If the estimation variance at distance h is less than or equal to 
the sample variance, then the estimated thickness at point (xth) using 
the thickness at point x is statistically reliable. Furthermore, the 
zone of influence determined from the variogram is also reliable. On 
the other hand, if the estimation variance is greater than the sample 
variance, the estimated thickness at point {xth) and zone of influence 
are unreliable. The following example will illustrate the situation. 

Suppose at point x, a borehole encounters a 5-foot thick seam. The 
geologist would like to know what is the thickness of this seam at, say, 
one mile away. From the data on the same seam in the area, he has cal- 
culated the sample variance as | foot (standard deviation = 1 fit) 

From the variogram, he obtains an estimation variance of 4 feet? (stan- 
dard deviation = 2 ft) at a distance of one mile. Assuming there is no 
sampling error, then the information from the variogram indicates that 

the seam thickness at a point one mile away from the borehole will be 
(542) feet. However, the actual range of the seam thickness at that point 
should be (541) feet as indicated by the sample variance. That means the 
variogram cannot be used to extrapolate seam thickness from one point to 
another with a distance of one mile in between. 

In the present study, the semi-variogram, v(h) shows a large nugget 

effect and crosses over the sample variance as the distance increases. 
As has already been discussed, the high nugget effect value decreases the 
reliability of the sample variance. Because of this unreliability, it is 
not possible to compare the sample variance with the estimation variance 
with any degree of certainty. Without this comparison, the zone (or 


range) of influence cannot be determined with confidence. Furthermore , 
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there is a possibility that the control point spacing in this study is 

not close enough to provide information on the initial portion of the 
variogram curve near the origin. It is this initial portion of the 
variogram that gives the vital evidence on the continuity of the deposit. 
Moreover, without geophysical logs of better resolution, the error intro- 
duced in the estimation of coal seam thickness rendered the statistics 
used in this variogram study unreliable. Therefore, closer well spacing 
as well as better log resolution are needed in order to determined whether 


the variogram can be used to predict palaeodrainage direction. 
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Chapter 9 
SUMMARY AND CONCLUSIONS 


The facies change from the distinct six-member Mannville Formation 
as defined by Nauss (1945) to the Undifferentiated Mannville Group (May- 
cock, 1967) occurs in the vicinity of Township 40 in west-central Sask- 
atchewan, and is attributed to the presence of the Battleford Arch. The 
circumstantial evidence for this is the pinching out of both the Dina and 
Cummings Members in that general vicinity and the absence of coal seam 8 
in Townships 37 to 39. The arch was high either during or after the 
deposition of seam 8 depending on whether the absence of the seam in Town- 
Ships 37 to 39 was due to non-deposition or erosion. Since no unconform- 
ity was detected on top of this seam, the absence of seam 8 in this area 
was probably due to non-deposition. 

The Battleford Arch also had a capale effect on the depositional 
environments of the Mannville Group. North of the arch, depositional 
environments changed from fluvial during Dina time to nearshore and/or 
deltaic in the later part of the Senne iit history. On the other hand, 
south of the Battleford Arch, Mannville sediments are mainly of contin- 
ental origin pee aieate by the lack of marine fossiis and the abundance 
of coal seams. 

Coal seams in the Mannville Group are scattered both geographically 
and stratigraphically as shown by their distribution maps, centre of 
gravity map and relative standard deviation map. Deposition of coal was 
~ confined mainly to the palaeotopographic lows during Mannville time 
(figures 8.1 and 6.4) with the loci of deposition migrating frequently, 


producing thin coal seams scattered vertically through the Mannville beds. 
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A total of ten seams were identified from geophysical borehole logs 
with only seams 3 and 8 having relative regional significance. Distri- 
bution of the minor seams can be explained by Price and Ball's (1971) 
model which postulates that the coal was of drift origin, perhaps forming 
from log jams on the river banks rather than originating in widespread 
Swampy areas. The two major seams, 3 and 8, may have been deposited in 
swamps on the fluvial plain, giving them larger areal distributions. 

The fact that seam 8 occurs both in the northern and southern portions of 
the study area makes it a useful stratigraphic marker for correlation 
within the Mannville Group. 

Coal resources in place calculated on a one-well-per-township basis 
were estimated to be 19,400 million (short) tons, including seams 2 feet 
thick or thicker. Using seam thickness as the variable, a variogram study 
was carried out in part of the area. Due to the lack of a reliable esti- 
mation of the sample variance for seam thickness, it was not possible to 
determine the zone (or range) of influence of a data point. With the 
knowledge of the zone of influence and the depositional environments of 
the study area, the variogram may be abe fu tool to predict palaeodrain- 
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